1 



o 
o 

(N 

in 



SDLCQ and String/Field Theory Correspondences* 

Uwe Trittmann, Department of Physics, Ohio State University, Columbus, OH, USA 

String/Field theory correspondences have been discussed heavily in recent years. Here, we describe a testing 
scenario involving a non-perturbative field theory calculation using the framework of supersymmetric discrete 
light-cone quantization (SDLCQ). We consider a Maldacena-type conjecture applied to the near horizon geometry 
of a Dl-brane in the supergravity approximation. Numerical results of a test of this conjecture are presented with 
orders of magnitude more states than we previously considered. These results support the Maldacena conjecture 
and are within 10-15% of the predicted results. We present a method for using a "flavor" symmetry to greatly 
reduce the size of the Fock basis and discuss a numerical method that we use which is particularly well suited 
for this type of matrix element calculation. Our results are still not sufficient to demonstrate convergence, and, 
therefore, cannot be considered to be a numerical proof of the conjecture. We update our continuous efforts to 
improve on these results and present some results on the way to higher dimensional scenarios. 
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1. Introduction 

The conjecture that certain field theories ad- 
mit concrete realizations as string theories on 
particular backgrounds has caused a lot of ex- 
citement in the last years. Originally, the so- 
called Maldacena conjecture 0| assured that the 
N = A supersymmetric Yang-Mills (SYM) the- 
ory in 3-1-1 dimensions is equivalent to Type IIB 
string theory on an AdS^ x background. Mean- 
while, other string/field theory correspondences 
have been conjectured. Attempts to rigorously 
test these conjectures have met with limited suc- 
cess, because our understanding of both sides of 
the correspondences is usually insufficient. The 
main obstacle is that at the point of correspon- 
dence, we require two conditions to hold which 
are mutually exclusive. Namely, we want a situa- 
tion where the curvature of the considered space- 
time is small in order to be able to use the super- 
gravity approximation to string theory. One also 
desires the corresponding field theory to be in a 
small coupling regime. So far it has been impossi- 
ble to find such a scenario. We present a way out 
of this dilemma by relaxing the second require- 
ment and performing a non-perturbative calcula- 
tion on the field theory side. To create a man- 
ageable situation, we chose a special string/field 
theory correspondence in order to apply the non- 



perturbative method, namely SDLCQ, at its op- 
timal working point. 

SDLCQ, or Supersymmetric Discretized Light- 
Cone Quantization, is a non-perturbative method 
for solving bound-state problems that has been 
shown to have excellent convergence properties, 
in particular in low dimensions. Therefore, we are 
looking for a (preferably) two-dimensional field 
theory, which is conjectured to be equivalent to 
a string theory. It turns out that the Yang-Mills 
theory with 16 supercharges in two dimensions 
1^ has its corresponding string theory partner in 
a system of Dl-branes in Type IIB string the- 
ory decoupling from gravity Q. Since both sys- 
tems have separately been studied in the litera- 
ture already, this systems is an optimal candidate 
to study the string/field theory correspondence. 

The next step is to find an observable that can 
be computed relatively easy on both sides of the 
correspondence. It turns out that the correla- 
tion function of a gauge invariant operator is a 
well-behaved object in this sense. We chose the 
stress-energy tensor T^" as this operator and will 
construct this observable in the supergravity ap- 
proximation to string theory and perform a non- 
perturbative SDLCQ calculation of this correla- 
tor on the field theory side. 



*Based on work with S. Pinsky, O. Lunin, and J.R. Hiller. 
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Figure 1. Phase diagram of two-dimensional M = (8, 8) SYM: the theory flows from a CFT in the UV to a 
conformal cr-model in the IR. The SUGRA approximation is vahd in the intermediate range of distances, 
^/gvAiVNc < X < ^/NllgvM- 



2. Correlation functions from supergravity 

It is instructive to take a closer look on the ex- 
pected properties oi M — (8, 8) SYM, before we 
proceed to technical details on the string theory 
side. In the extreme ultra-violet (UV) this the- 
ory is confornially free and has a central charge 
cuY — N'^ . Perturbation theory in turn will be 
valid for small effective couplings g — gyAiVT^x, 
where x is a space coordinate. For large dis- 
tances, in the far infra-red (IR), the theory be- 
comes a conformal a-model with target space 
{R^)^'^/Sn,. The central charge is cir = N^. It 
is a bit more involved to show that here pertur- 
bation theory breaks down when x^\/Nc/ gvM, 
see e.g. Ref. 

The intermediate region, 1/ gvAiVNc < x < 
V^c/gYM, where no perturbative field theoreti- 
cal description is possible, is fortunately exactly 
the region which is accessible to string theory; 
or rather, to the supergravity (SUGRA) approx- 
imation to Type IIB string theory on a special 
background. It is that of the near horizon geom- 
etry of a Dl-branc in the string frame, which has 
the metric 
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where we defined gvM = SP^^gYAiV^- In the 
description of the computation of the two-point 
function we follow Ref. ||]. The correlator has 
been derived in Ref. , being itself a generaliza- 
tion of Refs. 



First, we need to know the action of the di- 
agonal fluctuations around this background to 
the quadratic order. We would like to use the 
analogue of Ref. [D for our background, Eq. (||), 
which is not (yet) available in the literature. How- 
ever, we can identify some diagonal fluctuating 
degrees of freedom by following the work on black 
hole absorption cross-sections pO| , pT| . One can 
show that the fluctuations parameterized like 

ds' = {l + fix',U)+gix°,U))gooidxy 

+ {l + 5f{x'',U)+g{x'',U))gn{dx')' 

+ (1 + fix", U) + 5(x", U)) guudU^ 

+ + fix", U) - ^gix°, U)^ gnndn'r 

= il + 3fix°,U)-gix",U))e'f'°, (2) 
satisfy the following equations of motion 

f'iu) = -^.rm + ^fiu) (3) 

7 79 PP 

9"{U) = -Lg'iU) + l±giU) + ^-f^giU). 

Without loss of generality we have assumed here 
that these fluctuation vary only along the x° di- 
rection of the world volume coordinates, and be- 
have like a plane wave. One can interpret a Dl- 
brane as a black hole in nine dimensions. The 
fields /([/) and giU) are then the minimal and 
the fixed scalars in this black hole geometry. In 
ten dimensions, however, we see that they are re- 
ally part of the gravitational fluctuation. Conse- 
quently, we expect that they are associated with 
the stress-energy tensor in the operator field cor- 
respondence of Refs. mm. In the case of the cor- 
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respondence between M=A SYM field theory and 
string tiieory on an AdS^ x background, the 
superconformal symmetry allows for the identifi- 
cation of operators and fields in short multiplets 
|T2| . In the present case of a Dl-brane, we do 
not have superconformal invariance and this tech- 
nique is not applicable. Actually, we expect all 
fields of the theory consistent with the symmetry 
of a given operator to mix. The large distance 
behavior should then be dominated by the contri- 
bution with the longest range. The field U) 
appears to be the one with the longest range since 
it is the lightest field. 

Eq. (p|) for /([/) can be solved explicitly 
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where ^3/2 (2:) is a modified Bessel function. If 
we take f{U) to be the analogue of the minimally 
coupled scalar, we can construct the flux factor 



J- = lim 
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(5) 



up to a numerical coefficient of order one which 
we have suppressed. We see that the leading non- 
analytic contribution in is due to the term. 
Fourier transforming the latter yields 



{o{x)om = 



N2 

gvMX^ 



(6) 



This is in line with the discussion at the begin- 
ning of this section. We expect to deviate from 
the trivial scaling behavior of the corre- 

lator at xi = '^1 gYM^fN~c and X2 = VNc/gvM- 
This yields the phase diagram in Fig. 0. It is 
interesting to note that the entire Nc hierarchy 
is consistent in the sense of Zamolodchikov's c- 
theorem, which assures that the central charges 
obey c{x) > c{y), whenever x < y 

3. The correlator from SDLCQ 

Discretized Light-Cone Quantization (DLCQ) 
preserves supersymmetry at every stage of the 



calculation if the supercharge rather than the 
Hamiltonian is diagonalized [Q. The framework 
of supersymmetric DLCQ (SDLCQ) allows one 
to use the advantages of light-cone quantization 
{e.g. a simpler vacuum) together with the excel- 
lent renormalization properties guaranteed by su- 
persymmetry. Using SDLCQ, we can reproduce 
the SUGRA scaling relation, Eq. (^), fix the nu- 
merical coefficient, and calculate the cross-over 
behavior at l/gvAiVT^ < x < ^/Nc/gyM- To ex- 
clude subtleties, nota bene issues of zero modes, 
we checked our results against the free fermion 
and the 't Hooft model and found consistent re- 
sults. 

The technique of (S)DLCQ was reviewed in 
Ref. ||l^, so we can be brief here. The basic 
idea of light-cone quantization is to parameter- 
ize space-time using light-cone coordinates 



V2 



(7) 



and to quantize the theory making a;"*" play the 
role of time. In the discrete light-cone approach, 
we require the momentum p_ = p+ along the 
x~ direction to take on discrete values in units 
of / K where p+ is the conserved total mo- 
mentum of the system. The integer K is the 
so-called harmonic resolution, and plays the role 
of a discretization parameter. One can think of 
this discretization as a consequence of compact- 
ifying the x~ coordinate on a circle with a pe- 
riod 2L — 2-kK/p'^ . The advantage of discretiz- 
ing on the light cone is the fact that the dimen- 
sion of the Hilbert space becomes finite. There- 
fore, the Hamiltonian is a finite-dimensional ma- 
trix, and its dynamics can be solved explicitly. In 
SDLCQ one makes the DLCQ approximation to 
the supercharges Q^. Surprisingly, also the dis- 
crete representations of Qi satisfy the supersym- 
metry algebra. Therefore SDLCQ enjoys the im- 
proved renormalization properties of supersym- 
metric theories. To recover the continuum result, 
K has to go to infinity. Incidentally, what one 
finds is that SDLCQ usually converges faster than 
the naive DLCQ towards the continuum limit. 

Let us now return to the problem at hand. 
We would like to compute a general expres- 
sion for the correlator of the form F{x~ ,x'^) — 



4 



{O{x-,x+)O{0,0)). In DLCQ one fixes the total 
momentum in the x~ direction, and it is natural 
to compute the Fourier transform and express it 
in a spectrally decomposed form 

1 



2L 



{0{P-,x+)0{~P-,Q)) 



2L 

x(n|O(-P_,0)|0) 



(8) 



The form of the correlation function in position 
space is then recovered by Fourier transforming 
with respect to P_ — Ktt/L. We can continue 
to Euclidean space by taking r = •\/2x+x~ to be 
real. The result for the correlator of the stress- 
energy tensor is 

2 / 1 \ 2 
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(9) 



where Mi is a mass eigenvalue and K^lx) is the 
modified Bessel function of order 4. Note that 
this quantity depends on the harmonic resolution 
K, but involves no other unphysical quantities. 
In particular, the expression is independent of the 
box length L. 

The momentum operator T++(a;) of two- 
dimensional J\f — 8 SYM is given by 



T++{x) 



tr [{d-X^f 
+ - {iu^'d-u'^ -i(d-u'')u°') 



,(10) 



with /, a = 1 . . . 8. X and u are the physical ad- 
joint scalars and fermions, respectively j2j. When 
discretized, these operators have the mode expan- 
sions 



^1,3 - 
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The matrix element (L/7r)(0|T++(ii:)|i) can be 
substituted directly to give an explicit expression 
for the two-point function. We see immediately 
that the correlator has the correct small-r behav- 
ior, for in that limit, it asymptotes to 



N^i2nb + nf) 
47r2r4 



1 - 



1 

K 



which we expect for the theory of nb{nf) free 
bosons (fermions) at large K. 

On the other hand, the contribution to the cor- 
relator from strictly massless states is given by 



— F{x-,x+) 



E 



6 



(13) 



-(0|^++(if)|^) 

TT 



It is important to notice that this l/r"* behav- 
ior at large r is not the one we are looking for at 
large r. First of all, we do not expect any massless 
physical bound state in this theory, and, addition- 
ally, it has the wrong Nc dependence. Relative to 
the 1/r* behavior at small r, the l/r** behavior at 
large r that we expect is down by a factor of 1/iVc. 
This behavior is suppressed because we are per- 
forming a large- A'c calculation. All we can hope 
is to see the transition from the 1 /r'* behavior at 
small r to the region where the correlator behaves 
like l/r^. 

4. Symmetries and Numerics 

In principle, we can now calculate the corre- 
lator numerically by evaluating Eq. (|^). How- 
ever, it turns out that even for very modest har- 
monic resolutions, we face a tremendous numer- 
ical task. At K — 2,3,4, the dimension of the 
associated Fock space is 256, 1632, and 29056, re- 
spectively. The usual procedure is to diagonalize 
the Hamiltonian P~ and then to evaluate the pro- 
jection of each eigenfunction on the fundamental 
state T~^'^{—K)\0). Since we are only interested 
in states which have nonzero value of such pro- 
jection, we are able to significantly reduced our 
numerical efforts. 

In the continuum limit, the result does not de- 
pend on which of the eight supercharges Q~ one 
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chooses. In DLCQ, however, the situation is a 
bit subtler: while the spectrum of (Q^)^ is the 
same for all a, the wave functions depend on the 
choice of supercharge This dependence is an 
artifact of the discretization and disappears in 
the continuum limit. What happens if we just 
pick one supercharge, say Q^? Since the state 
T++(— iir)|0) is a singlet under R-symmetry act- 
ing on the "flavor" index of Q~ , the correlator (||) 
does not depend on the choice of a even at finite 
resolution! 

We can exploit this fact to simplify our cal- 
culations. Consider an operator S commuting 
with both and T++{-K), and such that 
S\Q) — sq\0). Then the Haniiltonian and S can 
be diagonalized simultaneously. We assume in the 
sequel that the set of states \i) is a result of such 
a diagonalization. In this case, only states satis- 
fying the condition S\i) = so\i) contribute to the 
sum in and we only need to diagonalize in 
this sector, which reduces the size of the problem 
immensely. We can deduce from the structure of 
the state T++(— _?C)|0) that any transformation of 
the form 



hUk)^g{a)by{k) 



(14) 



given arbitrary permutations P and Q of the 8 
flavor indices, commutes with T'^~^{—K). The 
vacuum will then be an eigenstate of this trans- 
formation with eigenvalue 1 . The requirement for 
P~ = (Qi)'^ to be invariant under S imposes 
some restrictions on the permutations. In fact, 
we will require that be invariant under S, in 
order to guarantee that P~ is invariant. 
The form of the supercharge from j^] is 



[...]bi{k3)aiik,)aiik2) + ... (15) 

+ (/3//?J - Pjl3j)a,f3[..]bl{k3)aiik,)aj{k2) + .... 

Here the /?/ are 8x8 real matrices satisfying 
= 26ij. We use the special represen- 
tation for these matrices given in Ref. . 

Let us consider the expression for Qj^, Eq. (p^. 
The first part of the supercharge does not include 
f3 matrices, and is therefore invariant under the 
transformation, Eq. ([T^), as long as g{l) = 1 and 



Q[l] — 1. We will consider only such transfor- 
mations. The crucial observation for the analysis 
of the symmetries of the (3 terms is that in the 
representation of the (3 matrices we have chosen, 
the expression Bfj — {PiPj — P,j(3j) , may take 
only the values ±2 or zero. Besides, for any pair 
(/, J) there is only one (or no) value of a corre- 
sponding to nonzero B. Using this information, 
we may represent B in a compact form. With the 
definition [§1 



r «, Bfj = 2 
flu ^ < ~a, Bfj = -2 

[ , Bfj = for all a 



(16) 



together with the special choice of (3 matrices 
we get the following expression for /i 
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The next step is to look for a subset of the 
transformations, Eq. (p^, which satisfy the con- 
ditions g{l) = 1 and Q[l] — 1 and leave the ma- 
trix /i invariant. This invariance implies that 



Q[Mp[/]p[.7]] .9(m/./)/(^)/(-^)m 



ij ■ 



(17) 



The subset of transformations we are looking for 
forms a subgroup R of the permutation group 
5*8 X S'g . Consequently, we will search for the ele- 
ments of R that square to one. Products of such 
elements generate the whole group in the case of 
5*8 X 5*8. We will show later that this remains 
true for R. Not all of the Z2 symmetries satisfying 
( p^ ) are independent. In particular, if a and b are 
two such symmetries then aba is also a valid Z2 
symmetry. By going systematically through the 
different possibilities, we have found that there 
are 7 independent Z2 symmetries in the group R. 
They are listed in Table 0. We explicitly con- 
structed all the symmetries of the type, Eq. (^4|), 
which satisfy Eq. (0) using Mathematica. It 
turns out that the group of such transformations 
has 168 elements, and we have shown that all of 
them can be generated from the seven Z2 sym- 
metries mentioned above. 
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Table 1 

Seven independent Z2 symmetries of the group R, which act on the 'flavor' quantum number of the 
different particles. Under the first of these symmetries, e.g., the boson oi is transformed into aj, etc. 



In our numerical algorithm we implemented the 
Z2 symmetries as follows. We can group the Fock 
states in classes and treat the whole class as a new 
state, because all states relevant for the correlator 
are singlets under the symmetry group R. as an 
example consider the simplest non-trivial singlet 

|l) = l^tr(at(l,/)at(X-l,/))|0). (18) 
1=1 

Hence, if we encounter the state a^ {l,l)a^ {K — 
1,1) |0) while constructing the basis, we will re- 
place it by the class representative; in this case, 
by the state |1). Such a procedure significantly 
decreases the size of the basis, while keeping all 
the information necessary for calculating the cor- 
relator. In summary, this use of the discrete fiavor 
symmetry of the problem reduces the size of the 
Fock space by orders of magnitude. 

In addition to these simplifications, one can fur- 
ther improve on the numerical efficiency by using 
Lanczos diagonalization techniques. Namely, we 
substitute the explicit diagonalization with an ef- 
ficient approximation. The idea is to use a sym- 
metry preserving (Lanczos) algorithm. If we start 
with a normalized vector \ui) proportional to the 
fundamental state T^~^{—K)\0), the Lanczos re- 
cursion will produce a tridiagonal representation 
of the Hamiltonian Hlc = 2P'^P~. Due to or- 
thogonality of only the (1,1) element of 
the tridiagonal matrix, i?i.i, will contribute to 
the correlator. We exponentiate by diagonalizing 
HlcVz = Kvi with eigenvalues Ai and get 



Finally, we Fourier transform to obtain 

Nl 

X E |(i'j)il'A,'i^4(v/2^+^), 

which is equivalent to Eq. (^). This algorithm is 
correct only if the number of Lanczos iterations 
Nl runs up to the rank of of original matrix. But 
in praxi already a basis of about 20 vectors covers 
all leading contribution to correlator [ p^ . 

5. Results 

To evaluate the expression for the correlator 
J'(r), we have to calculate the mass spectrum and 
insert it into Eq. (^. In the M = (8,8) super- 
symmetric Yang-Mills theory the contribution of 
massless states becomes a problem. These states 
exist in the SDLCQ calculation, but are unphys- 
ical. It has be shown that theses states are not 
normalizable and that the number of partons in 
these states is even (odd) for K even (odd) ||^. 
Because the correlator is only sensitive to two par- 
ticle contributions, the curves J^{r) are different 
for even and odd K. Unfortunately, the unphys- 
ical states yield also the typical 1/r* behavior, 
but have a wrong Nc dependence. The regular 
l/r^ contribution is down by 1/Nc, so we cannot 
see this contribution at large r, because we are 
working in the large Nc limit. 

We can use this information about the unphys- 
ical states, however, to determine when our ap- 
proximation breaks down. It is the region where 
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Figure 2. Top: (a) Log-Log plot of T[r) = 
(r++(x)T++(0))(f^)',,^^^ vs. r for 

dYM^cf-n- = 1.0, K = 3,4,5 and 6. Bottom: 
(b) the log-log derivative with respect to r of the 
correlation function in (a). 



the unphysical massless states dominate the cor- 
relator sum. Unfortunately, this is also the re- 
gion where we expect the true large-r behavior to 
dominate the correlator, if only the extra states 
were absent. In Fig. ||(a) for even resolution, the 
region where the correlator starts to behave like 
l/r"* at large r is clearly visible. In Fig. ^(b) we 
see that for even resolution the effect of the mass- 
less state on the derivative is felt at smaller values 
of r where the even resolution curves start to turn 
up. Another estimate of where this approxima- 
tion breaks down, that gives consistent values, is 
the set of points where the even and odd reso- 
lution derivative curves cross. We do not expect 
these curves to cross on general grounds, based 
on work in , where we considered a number of 
other theories. Our calculation is consistent in 
the sense that this breakdown occurs at larger 
and larger r as K grows. 

We expect to approach the line dT{r)/dr = — 1 
line signaling the cross-over from the trivial 
behavior to the characteristic behavior of 

the supergravity correlator, Eq. (||). Indeed, the 
derivative curves in Fig. |^(b) are approaching —1 
as we increase the resolution and appear to be 
about 85 — 90% of this value before the approx- 
imation breaks down. There is, however, no in- 
dication of convergence yet; therefore, we cannot 
claim a numerical proof of the Maldacena con- 
jecture. A safe signature of equivalence of the 
field and string theories would be if the derivative 
curve would flatten out at —1 before the approx- 
imation breaks down. 

6. Conclusions 

In this note we reported on progress in an at- 
tempt to rigorously test the conjectured equiva- 
lence of two-dimensional N — (8, 8) supersym- 
metric Yang-Mills theory and a system of Dl- 
branes in string theory. Within a well-defined 
non-perturbative calculation, we obtained results 
that are within 10-15% of results expected from 
the Maldacena conjecture. The results are still 
not conclusive, but they definitely point in right 
direction. Compared to previous work we in- 
cluded orders of magnitude more states in our 
calculation and thus greatly improved the test- 
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Figure 3. Log-Log plot of the three- 
dimensional correlation function / = 

r^T++{x)T++{0)) (s^Y ^^^^ vs. r for 

g = g'^j^.jy/NJ/2TT^/^ = 1.0 for K = 6 and T = 1 
to 5. 

ing conditions. We remark that improvements of 
the code and the numerical method are possible 
and under way. During the calculation we no- 
ticed that contributions to the correlator come 
from only a small number of terms. An analytic 
understanding of this phenomenon would greatly 
accelerate calculations. We point out that in prin- 
ciple we could study the proper l/r behavior at 
large r by computing 1/Nc corrections, but this 
interesting calculation would mean a huge numer- 
ical effort. 

7. Outlook 

It remains a challenge to rigorously test the 
conjectured string/field theory correspondences. 
Although the so-called Maldacena conjecture 
maybe the most exciting one, because it promises 
insight into full four-dimensional Yang-Mills theo- 
ries in the strong coupling regime, there are other 
interesting scenarios. For instance, it was conjec- 
tured that the supergravity solutions correspond- 



ing to p -I- 1 SYM theories are black p-brane so- 
lutions, see e.g. Ref. ||]. Consequently, there 
are interesting testing scenarios also in three- 
dimensional spacetime. Numerically, of course, 
things get the more difficult, the more dimen- 
sions are involved. On the way to the full four- 
dimensional problem, it may be worthwhile to 
present our latest results on correlation functions 
in three dimensions, see also Fig. |^ shows 

the correlator for M —1 SYM(2-fl) as a function 
of the distance r: it is converging well with the 
transverse cut-off T . To put things in perspec- 
tive, we note that the largest Hamiltonian matrix 
involved in this calculations requires to set up a 
Fock basis of approximately two million states. 
This is by a factor 100 more than we used in the 
test of the Maldacena conjecture described in this 
article, which itself was already substantially bet- 
ter than the first feasibility study ||^. 

We hope to proceed on this way and to be able 
to present interesting results soon. 
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